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While both the hole and electron doped cuprates can exhibit d x 2 _ y i- wave super¬ 
conductivity, the local distribution of the doped carriers is known to be significantly 
different with the doped holes going primarily on the O sites while the doped electrons 
go on the Cu sites. Here we report the results of density-matrix-renormalization- 
group calculations for a three-orbital model of a Cu02 lattice. In addition to the 
asymmetric dependence of the intra-unit-cell occupation of the Cu and O for hole and 
electron doping, we find important differences in the longer range spin and charge 
correlations. As expected, the pair-field response has a d x 2 _ y 2 -\ike structure for both 
the hole and electron doped systems. 

How well does a 3-orbital Hubbard model describe the properties of the cuprates? These 
materials are known to be charge-transfer systems and from the analysis of Zaanan, Sawatzky 
and Allen 1 -' one would expect that a minimal model which includes a Cu 3d x 2_ y 2 orbital and 
two O 2pcr orbitals per unit cell would be required. Indeed, early on a 3-orbital Hubbard 
model was proposed by several groups 2 ^, and various quantum Monte CarlcP 5! 'and embedded 
cluster calculation^ 6 7 have shown that this model exhibits a number of the basic magnetic 
and single particle spectral weight properties that are seen in the cuprates. More recently, 
experimental measurements of both the hole and electron doped cuprates have provided new 
information on the spatial charge and spin structure which can occur when these materials 
are dopecP®. So the question of whether a 3-orbital Hubbard model provides a suitable 
framework with which to describe the physics of the cuprates has been enlarged. Here 
with the experimental results for hole doped La 2 - x Ba x Cu 04 (LBCO) and electron doped 
Nd 2 - a; Ce a ;Cu 04 (NCCO) in mind, we have carried out density matrix renormalization group 
(DMRG) 1 -^ calculations with the goal of determining whether the 3-orbital Hubbard model 
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remains an adequate model for the cuprates. 

Neutron scattering studies of the LTT phase of LBCO find that the doped holes form a 
striped structure consisting of regions with excess holes separated by 7r-phase shifted anti¬ 
ferromagnetic region#®. At 1/8 hole doping, superconducting correlations are observed to 
onset together with the stripe order 11 . This pair-density-wave phase is believed to have a 
d -wave pair-field which is large in the regions with excess holes and oscillates in sign be¬ 
tween these charged regions 12 13 . Achkar et al. 14 have reported resonant soft x-ray scattering 
measurements which show that the charge distribution on the oxygens of LBCO have an 
s'-CDW orbital structure in which the charge modulations on the 0 Px and 0 Py sites in a 
unit cell are in phase. STM Studies of BSCCO (p ~ 8%) and NaCCOC (p ~ 12%) find 
that these materials have a predominantly d-CDW orbital form factor in which these 0 Px 
and 0 Py charge modulations are out of phas#P Finally, recent resonant x-ray scattering 
measurements of Nd^zCe^CuCR near optimal doping 16 find charge order which occurs with 
a similar periodicity and Cu-0 bond orientation to that of the charge stripes seen in LBCO. 

One-band Hubbard and t-J models have been found, within various approximations, to 
exhibit striped charge and spin structure#®®} modulated nematic phase#®"® a s well as 
pair density wave phased® 30 . RPA calculations for the three-band Hubbard model have 
also found nematic phases in certain parameter regimes 31 33 . Earlier DMRG calculations for 
a 3-orbital model of a two-leg Cu0 2 ladder showed the expected local asymmetric charge- 
transfer behavior in which doped holes tend to predominantly go on the 2p<r orbitals while 
doped electrons go on the Cu 3d x 2_ y 2 orbital# 4 35 l These calculations also found d x 2 _ y 2 -like 
pairing correlations for both hole and electron doping in which the near neighbor Cu rung 
and leg pair-field correlations differ in sign. Here we extend these calculations to an 8 x 4 
Cu0 2 cluster with cylindrical boundary conditions. The cylindrical boundaries reduce the 
edge effects associated with the ladder, more reliably representing bulk behavior. The Lx4 
geometry is also the minimal size that can contain stripe-like clusters of holes. With the 
8x4 system we study the tendencies towards striping in the hole densities and whether 
doped holes or electrons modulate the phase of the antiferromagnetism. We also study the 
hopping kinetic energy associated with added holes or electrons, and the pairing tendencies 
in the doped system. 

The lattice structure and the parameters of the three orbital Cu0 2 model that we will 
study are shown in Fig. [lj The model has a Cu0 2 unit cell consisting of a 3d x 2_ y 2 orbital on 
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FIG. 1: The C 11 O 2 lattice with onsite Cu ?>d x i_ y i and O 2 p x and 2 p y energies and e p , near 
neighbor Cu-0 and 0-0 hoppings t p d and t pp , onsite Cu and O Coulomb interactions Ud and U p and 
near neighbor Cu-0 and 0-0 Coulomb interactions V p d and V pp , respectively. The Hamiltonian 
for this 3-orbital model is given in a hole representation by Eq. ([!]). Here we will work with energies 
measured in units of t p d, and take as a generic set of parameters t pp = 0.5, A p d = e p — = 3, 

U d = 8, U p = 3, V pd = 1 and V pp = 0.75. 

the Cu site and 2p x /2p y orbitals on the x and y oxygens. I 11 a representation in which the 
vacuum of the 3-orbital model has the configuration ( d 2 x2 _ y 2 p 2 x pl ), d J. and p+ a create holes 
with spin a on the i th Cu and j th O sites respectively, and the Hamiltonian has the form 

H = A p d "Y, Pin-Pier ~ tpd 'Y, (d ia pj a + Pj a di a ) 
ia ( ij)<r 

~ Y. iptrPjo- + PjaPia) 

+ u d Y n % + u pY n A n \ 

+ Vpd rfrij + Vpp nWj (!) 
(u> <b> 

Here A ps = e p — is the energy difference between having a hole on an O site versus a Cu 
site, t p d and t pp are one-hole hopping matrix elements between near-neighbor Cu and O sites 
and near-neighbor O sites, respectively. The sums (ij) in Eq. [ljdenote sums over the relevant 
nearest-neighbor sites. Ud and U p are the onsite Cu and O Coulomb interactions and V p d 
and Vp P are the nearest-neighbor Cu-0 and 0-0 Coulomb interactions, respectively. The 
phases of the orbitals have been fixed such that the signs of the hopping matrix elements 
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remain the same throughout the lattice and are positive. 

The hopping parameters t pd and t pp found for La 2 Cu0 4 and Nb 2 Cu0 4 in various cluster 
and LDA calculations are relatively close to each other. We will measure energies in units 
of t pdj and set t pp /t pd = 0.5 for both of these materials® 3 ® 1 . The primary difference in the 
one-electron parameters occurs in A pd where the absence of the apical oxygens in Nd 2 Cu0 4 
is expected to lead to a reduction in A pd relative to La 2 Cu0 4 . Indeed this is found in LDA 
calculations, however the appropriate bare values of A pd to use in the 3-orbital Hamiltonian 
Eq. ([!]) has posed a problem because of double counting corrections 63 '. Here we find that 
setting A pd /t p d = 3 gives reasonable values for the charge gap and exchange interaction. 
Thus, working in units of t p d we will take for a canonical set of parameters 

t pp = 0.5, A pd = e p - e d = 3, U d = 8, Up = 3, V pd = 1, V pp = 0.75 (2) 

These parameters are appropriate for a charge transfer system for which U d > e p — £ d = A pd 
and Apd > 2 t pd . Using a similar set of parameters for a 2-leg Cu0 2 ladder we previously 
found at half-filling a charge gap A c ~ t pd and a spin gap A s ~ 0.03f p d- For a 2-leg ladder the 
effective exchange coupling J ~ 2A S ~ 0.06t pd . For t pd of order 1 to 2 eV, these correspond 
to reasonable values for the charge gap and the exchange interaction. Our plan is to use 
this same set of parameters for both the hole and electron doped systems and focus on the 
differences that arise between them. We will comment on the effect of reducing t pp and, for 
the electron doped case, the effect of reducing A pd . 

The DMRG calculations will be carried out for an 8 x 4 Cu0 2 lattice which has periodic 
boundary conditions in the 4-unit cell y-direction and open ends in the 8-unit cell x-direct ion. 
For the charge and spin studies, we will work with a fixed number of holes 32 + A and a 
hole density per Cu0 2 unit cell x = 1 + A/32 which is 1 for the undoped system. Positive 
values of A (x > 1) correspond to hole doping and negative values of A (x < 1) to electron 
doping. We typically did 15 DMRG sweeps, keeping up to m = 4000 states on the last 
sweep. This led to excellent convergence for the local quantities that we report here. A 
typical maximum truncation error was ~ 10 -5 ; extrapolating the truncation error to zero 
gave typical fractional errors in the total energy also about ~ ICR 5 . Without extrapolation, 
fractional errors in energy were estimated to be less than ICR 4 , and absolute errors in local 
quantities were in the range ICR 3 — ICR 4 . The good overall convergence for this cluster 
suggests that wider systems, say up to width 6, will be accessible for near-future studies. 
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In Fig. [2] we show the effect of doping on the local charge density and squared spin 
moments on the Cu and O sites as a function of the hole density x. As we will discuss 
later, there can be inter- and intra-cell spatial structure in the charge and spin. The results 
shown in Fig. [2] represent site averages taken over the 8x4 lattice. For the undoped x — 1 
(N = 0) system where there is one hole per C 11 O 2 unit cell, Fig. [2] shows that the hole 
occupation is approximately 80% on the Cu site and 10% on each of the two O sites for the 
parameters that we have chosen. When additional holes are added they go approximately 
75% onto the two O sites and 25% onto the Cu site of the unit cell. Alternatively, under 
electron doping, the added electrons go approximately 90% onto the Cu site and only 10% 
onto the two O sites. This is of course what one would expect for a charge-transfer system. 
The change in the square of the spin moments on the Cu and O sites is seen to vary with 
the hole concentration x in a similar manner to that of the charge occupation. For electron 
doping (x < 1), an electron added to a Cu site removes the hole spin moment leading to a 
decrease in (S' 2 ) averaged over the lattice, while for hole doping the square of the O hole 
spin moment increases as holes are primarily added to the O sites. 

To study the longer range spin and charge correlations, we have applied a weak staggered 
magnetic field to the Cu sites on the left hand edge of the 8x4 lattice. The expected 
antiferromagnetic response of the undoped system is shown in Fig. [3j Here, the diameters 
of the circles are proportional to the density of the holes and the lengths of the arrows are 
proportional to the spin moments. One sees, as shown in Fig. [2] that the holes are mainly 
on the Cu sites. The applied edge held has broken the spin symmetry and there is a well 
formed antiferromagnetic spin pattern. 

In Fig. [4] we contrast the results for hole doping on the left with electron doping on the 
right. In this figure, the hole density distribution of the undoped lattice shown in Fig. [3] has 
been subtracted. The diameters of the circles for hole doping on the left are proportional 
to the added hole density while the diameters of the circles on the right are proportional 
to the added electron density. In this figure the diameter scale used for the hole density is 
0.1 and for the electron density 0.15. In the top left hand lattice shown in Fig. [4| 2 holes 
have been added to the 32 holes of the undoped 8x4 lattice giving x = 1.0625. The lattices 
shown below this have 4,6 and 8 holes added corresponding to hole concentrations x per 
Cu0 2 unit cell of 1.125, 1.1875 and 1.25, respectively. The lattices on the right hand side 
of Fig. [I] show similar results for the case in which electrons are added (or holes removed). 
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FIG. 2: The hole density (n) and the square of the spin moment ( S 2 ) = | (rq- + n i — 2n^n^) on 
the Cu (black circles) and O (red squares) sites versus the hole density x = 1 + A/32 per Cu02 
unit cell. The undoped 8x4 CuC >2 lattice has 32 holes and x = 1. N > 0 (x > 1) corresponds 
to doping additional holes while N < 0 (x < 1) corresponds to electron doping which reduces the 
number of holes. 

From top to bottom these lattices have 30, 28, 26 and 24 holes, respectively, corresponding 
to x = 0.9375, 0.875, 0.8125 and 0.75. As in Fig. [3j a staggered magnetic field (h = 0.1) 
was applied to the Cn sites on the left-hand edge of the lattice. 

As shown in Fig. [2j the additional holes tend on average to go onto the O sites, but as 
seen in Fig. [4] their distribution is not uniform. For hole doping there is a tendency for 
stripe formation separated by 7r-phase shifts in the antiferromagnetic correlations. For x = 
1.125 there are two, approximately Cu-site centered, stripes separated by a 7r-phase shifted 
antiferromagnetic region similar to the well known behavior of Lai^gNdo^Sro.^C-viCXi 10 - For 
x = 1.0625 there is a single stripe and for x = 1.1875, corresponding to the addition of 6 
holes on the 8x4 lattice, one can see the remnants of a three stripe structure. This structure 
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FIG. 3: The hole occupation (n) and spin structure (S z ) for the undoped (32 hole) 8x4 Cu02 
lattice. The hole occupation is proportional to the diameter of the circles. A staggered magnetic 
field of magnitude h = 0.1 was applied to the Cu sites along the left-hand edge of the 8x4 Cu02 
lattice which has periodic boundary conditions in the y direction and open end boundary conditions 
in the x-direction. 

vanishes for the strongly overdoped x = 1.25 case. The stripe spacing for the three lower 
hole dopings is consistent with the relation dr 1 = 2(x — 1) and the well known spin 5 spin and 
charge Charge incommensurability relation Charge = 2<5 spin found in the La-based cuprate^®. 
A closer look at the structure of the charge and spin distributions for the 1/8 (x — 1.125) hole 
doped lattice is shown in Fig. [5j)a) . Here a weak staggered magnetic held has been applied 
to both ends of the 8x4 lattice. In this case, the x = 1.125 hole doped system exhibits 
bond centered charged stripes separated by 7r-phase shifted antiferromagnetic regions. The 
charge modulations on the 0 Px and 0 Py sites are in phase leading to what was called an 
s'-CDW-SDW phase in Ref. [Til There may be an additional small admixture of d-CDW. 
Of course the 8x4 lattice already breaks C4 symmetry so one expects differences in the x 
and y oxygen hole occupations. Increasing V pp leads to an increase in these differences 38 but 
the s' symmetry remains dominant. 

For the electron doped system, one sees on the right hand side of Fig. [4] that the spin 
and charge structure appears quite different from the hole doped case. Of course the doped 
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FIG. 4: Charge and spin structure for the hole (left) and electron (right) doped lattices for the 
parameters given in Eq. <©• Here the hole density distribution for the undoped lattice shown in 
Fig. | has been subtracted. The added hole density on a site is proportional to the diameter of the 
red circles shown on the left. Similarly, the diameter of the blue circles on the right is proportional 
to the added electron density. Note the difference in the diameter scales for the hole and electron 
doped figures. The small ~ 5% of the additional electron density that goes onto an O site is not 
visible on this scale. For the 8x4 lattice, the left hand figures (top to bottom) correspond to the 
addition of 2,4,6 and 8 holes respectively, while the right hand figures correspond to the addition 
of a corresponding number of electrons. A staggered magnetic field of magnitude h = 0.1 has been 
applied to the Cu sites on the left-hand edge of the lattice. 
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electrons go dominantly on to the Cu sites and initially the small concentration of added elec¬ 
trons are repelled from the open edge boundaries by the “infinite” edge potential. However, 
for these parameters, by the time the electron doping reaches 0.125, a relatively uniform 
density of the added electrons is spread over the Cu sites and the antiferromagnetic or¬ 
der remains. A closer look at the 1/8 electron doped lattice is shown in Fig. |5|b). Here 
one can see that there are two charge stripes but the antiferromagnetic correlations remain 
commensurate. Thus for these parameters we find charge stripes with incommensurate an¬ 
tiferromagnetism for hole doping and commensurate antiferromagnetic spin correlations for 
electron doping. This remains the case for the electron doped system when Ais reduced 
as is expected in the T' structure where the apical oxygens are absent. Another important 
parameter is t pp which determines the effective hopping t' between next-near neighbor Cu 
sites. In Hubbard and t — t! — J models it is known that t' affects the stripe stability 1221 . 
Here we find that when the oxygen-oxygen hopping t pp is reduced, the amplitude of the 
charge stripes is increased and the spin structure for the electron doped system also be¬ 
comes incommensurate as shown in Fig. |5](c) for t pp = 0. The effect of reducing t pp acts to 
increase the frustration associated with the antiferromagnetic background and gives rise to 
the 7r-phase shifted antiferromagnetic regions separating the charge stripes. We find that 
when t pp is reduced (below < 0.25), striping can occur for both the electron and hole doped 
system. However, the tendency for striping is stronger in the hole doped system. 

The addition of holes to the filled band vacuum configuration ( d 2 x 2 _ y 2 p 2 x p ^) lowers the 
kinetic energy. For the hole Hamiltonian, Eq. ([!]), with positive hopping parameters this 
means that one expects the Cu-0 and 0-0 hole hopping strengths J2 s {djsPj S + p^ s di s ) and 
EsipfsPjs + P^sPis) to be positive for the 32 hole doped system as illustrated in Fig. J6|a). 

The Cu-0 hopping strength is larger than the 0-0 hopping strength reflecting the fact 
that the doped holes are of order 80% on the Cu sites. When additional holes are added, 
the hopping strength increases further. The difference in the hopping strengths between the 
36 hole doped lattice and the undoped 32 hole lattice are illustrated in Fig. |6](b). In the 
case of hole doping, the holes are distributed to both the Cu and O sites (~ 25% to the Cu 
and ~ 37.5% to each of the O sites) leading to the enhancement of both the Cu-0 and 0-0 
hopping strengths shown in Fig. |6](b). In addition, one sees evidence of the charge stripe 
structure. 

For the case of electron doping, the electrons go dominantly on the Cu sites. This 
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FIG. 5: (a) The charge and spin structure of the 1/8 (x = 1.125) hole doped system with a weak 
staggered magnetic held h = 0.1 applied to both ends of the 8x4 lattice. Here an s'-CDW- 
SDW structure is seen, (b) The charge and spin structure of the 1/8 (x = 0.875) electron doped 
system with t pp = 0.5. Here we find only a weak charge modulation and a commensurate spin 
antiferromagnetic structure, (c) Similar to the x = 0.875 electron doped system shown in (b) but 
with t pp = 0.0. In this case there is an incommensurate antiferromagnetic spin structure similar to 
that of the hole doped system. 

reduction of the average number of holes on the Cu sites leads to the reduction in the Cu-0 
hole hopping strength as shown in Fig. |6|c). Although the reduction of the average hole 
occupation on each O is only of order 5%, one might have expected that this would also 
reduce the strength of the 0-0 hole hopping. However, as shown in Fig. [6^c), the 0-0 
hopping strength is in fact slightly increased. The overall change in hopping strength is 
significantly smaller for the electron doped system. The total change in the kinetic energy 
measured in units of t p d per added hole is of order —3.2 while per added electron it is only 
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FIG. 6: (a) The hopping strength on each bond for the undoped system, (b) The difference in the 
hopping strength relative to the undoped system for the hole doped system with x = 1.125. (c) 
The same as (b), but for the electron doped system, x = 0.875 
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+0.7. This is consistent with the notion that the doped holes will enter a region of the band 
between T and M where there is significant dispersion while the electrons will enter near X 
where the dispersion is flat. 

In order to study the pairing response, we have applied a proximity singlet pair-held that 
couples to near neighbor Cu sites along the x direction, 

E (a+(4,4) + M4,4)) (3) 

(4 M 

with A x (4,4) = (d4 + i,4t rf 4,4+ - ^4+1,41^4,4t) /v^- The Cu-0 near neighbor responses 
in the ^-direction (A x (£ x , i y ) + A+ (4,4))/2 has a negative sign and is shown as the dashed 
lines in Fig. [7} The solid lines, which indicate a positive value, show the pair-held response 
(A„(4,4) + A+(4,4))/2 between near neighbor Cu-0 sites in the y direction. For these 
calculations the average hole number was set by a chemical potential //. In the top panel, 
Fig. 7(a), /i = —1.0, giving an average hole number (N) = 36.05 (x ~ 1.125) while for the 
lower panel, Fig. 7(b), /i = —3.2 giving (N) = 28.14 (x ~ 0.875). Both the hole doped 
pair-held response shown in Fig. 7(a) and the electron doped case shown in Fig. 7(b) have 
the expected d-wave-like sign change. 

If the proximity pair-held is applied only between the horizontal Cu-Cu sites on the left 
edge of the lattice, the induced pair-held decays rapidly in the x-direction for the hole doped 
system and somewhat more slowly for the electron doped case. The longer range pair-held 
correlations are suppressed by finite size effects. These are particularly severe for the periodic 
in y (tube-like) geometry of our Cu0 2 lattice. As seen in Fig. [5] for the hole doped lattice, 
a stripe appears each time a pair of holes is added for 2, 4 and 6 holes. This is consistent 
with previous 2-leg ladder studies where it was found that the preferred filling was 2 holes 
per 4 rungf 23 21 . Thus for an 8 x 4 Cu0 2 tube a low-energy fluctuation of ±2 holes involves 
the creation or destruction of a stripe, leading to a high energy spin configuration with a 
domain wall without holes. Alternatively, one could consider a configuration which has 4 
holes in a stripe, but this is also energetically unfavorable. While this effect is less severe for 
the electron doped system shown in Fig. 7(b) and the pair-field response is stronger because 
it lacks the antiferromagnetic domain walls, we expect that the small 8x4 size of the lattice 
still acts to suppress the hole number fluctuations. Properly comparing the pairing between 
the electron doped and hole doped systems will require larger systems, and if there are 
stripes they should be long, running lengthwise down the cylinder, either horizontally or 
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FIG. 7: The pair-field induced by the proximity pair-field, Eq. (3), for the hole x ~ 1.125 (a) 
and electron x ~ 0.875 (b) doped systems. Here the dashed lines denote a negative pair-field 
amplitude between near neighbor Cu-0 pairs of sites while the solid lines denote a positive pair- 
field amplitude. The applied proximity pair-field is between near neighbor Cu-Cu pairs of sites in 
the x-direction with magnitude Ao = 0.2, and the response shows the expected ( 4,2 2 -like behavior 

for both the hole and electron doped systems. 

spiraling. 

In summary, we have studied an 8 x 4 three orbital Hubbard model for C 11 O 2 with pa¬ 
rameters chosen to give a realistic charge gap and exchange coupling. With one hole per 
Cu0 2 unit, the hole occupation is approximately 80% on the Cu d x 2_ y 2 orbital and in the 
presence of a weak staggered edge magnetic field commensurate antiferromagnetic corre¬ 
lations are found to extend across the lattice. When additional holes are added they go 
~ 75% onto the O sites and charge stripes separated by 7r-phase shifted antiferromagnetic 
regions appear. The O hole occupation and the Cu spin structure has an s'-CDW-SDW 
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like structure 1 11 . When additional electrons are added, they go approximately 90% onto 
the Cu d x 2_ y 2 orbitals. For our small cluster, there is a weak tendency for charge modu¬ 
lations but the antiferromagnetic spin correlations remain commensurate. However, when 
the oxygen-oxygen one electron hopping t pp is reduced, a clear striped structure appears 
with incommensurate antiferromagnetic correlations. For both the hole and electron doped 
systems, the response of the y-near-neighbor Cu-Cu pair-held is out of phase (cl -wave like) 
with respect to the .x-near neighbor Cu-Cu pair-held induced by an applied x-near-neighbor 
proximity pair-held. These pair-held correlations are short range reflecting the finite size 
and geometric restrictions of the Cu0 2 cluster studied. 
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